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We study the role of Penning ionization on the photoassociation spectra of He(2'^Si)-E[e(2'^Si). 
The experimental setup is discussed and experimental results for different intensities of the probe 
laser are shown. For modelling the experimental results we consider coupled-channel calculations of 
the crossing of the ground state with the excited state at the Condon point. The coupled-channel 
calculations are first applied to model systems, where we consider two coupled channels without 
ionization, two coupled channels with ionization, and three coupled channels, for which only one of 
the excited states is ionizing. Finally, coupled-channel calculations are applied to photoassociation 
of He(2"^Si)-He(2"^Si) and good agreement is obtained between the model and the experimental 
results. 

PACS numbers: 33.20.-t,33.70. Jg, 34.50.Gb 
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I. INTRODUCTION 

Photoassociation spectroscopy (PAS) is a powerful 
technique for the study of ultracold collisions. In the 
case of alkali-metal atoms, the technique has been used to 
study long-range molecular states, the scattering length 
of the interacting atoms, the lifetime of the resonant 
state, and the dynamics of the interacting system [3| . In 
PAS two interacting ground-state atoms absorb a photon 
from the laser beam and are excited to an excited state 
of the molecule. The detection of the PAS process is in 
many cases performed by observation of trap loss, where 
one of the interacting atoms gains sufficient energy in the 
process to be lost from the trap. Since the relative energy 
at low temperatures is small, PAS is able to produce very 
high resolution spectra (<1 MHz). This allows for the de- 
tailed study of singly excited molecular states, where the 
interaction takes place at long range close to the dissoci- 
ation limit. Since the interactions at these distance are 
caused by dispersion forces, which are determined mainly 
by atomic parameters, detailed comparison between the- 
ory and experiment is possible for such systems. Since 
the strength of the PAS resonances depends on the over- 
lap between the ground and the singly excited state at 
the Condon point, the progression of the intensities of 
PAS lines allows for a determination of the wavefunc- 
tion of the ground state and thus for a determination of 
the scattering length, which is important for the achieve- 
ment of Bose-Einstein condensation. Due to this feature 
PAS has become a proven technique for interacting alkali- 
metal atoms. 

In the case of metastable rare-gas atoms, the situa- 
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tion is different. Two interacting, metastable rare-gas 
atoms have sufficient energy for the ionization of one 
of the atoms leading to the well-known process of Pen- 
ning ionization. The excitation of singly excited states 
in such systems can be accompanied by a strong increase 
of Penning ionization, which may prohibit the observa- 
tion of narrow resonances. Although the detection of ions 
and/or electrons can be much more sensitive than the de- 
tection of fluorescense as a measure of trap loss, observa- 
tion of PAS in metastable rare-gas systems using ion de- 
tection might not be the appropriate method. However, 
PAS has been observed a couple of years ago using ion 
detection for interacting, metastable helium atoms 
Since that time studies have been carried out for this sys- 
tem using both ion and fluorescense detection ^ . 
Recently, the results of two groups for PAS on metastable 
helium atoms have been compared, where one group used 
ion detection and the other group applied trap loss de- 
tection [^. The study showed, that the results of these 
two different techniques agreed within experimental reso- 
lution. This is especially remarkable, since both the tem- 
perature and the densities between the two experiments 
differed by three orders of magnitude. 

In this paper we study the role of Penning ionization on 
the PAS process. After discussing the technique of PAS 
for metastable helium and the results of our experiment, 
we introduce a pedagogical model for the PAS process, 
where the ground and excited channels are coupled by the 
light interaction. In this model we first consider the inter- 
action of one ground-state channel with one excited state 
channel without ionization. We show how the light inter- 
action leads to the observation of the resonances in the 
elastic scattering rate, which can be ascribed to bound 
states in the excited-state potential. We also show that 
by increasing the Franck-Condon overlap between the 
ground and the excited state, the resonances broaden and 
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shift due to the interaction with hght. If we allow the ex- 
cited state to become ionizing, the bound states obtain 
a finite lifetime and now also lead to resonances in the 
ionization cross section. The results of this two-channel 
problem are compared to a semi-classical Landau-Zener 
model. We find good agreement, if we optimize the pa- 
rameters of the Landau-Zener model with respect to the 
coupled-channel calculations. In a next step we expand 
our model by including a second, excited state and as- 
suming that only one of the excited states is ionizing. 
This leads to the observation of PAS resonances of the 
non-ionizing channel through the coupling with the ion- 
izing channel. Again, the results are compared with a 
Landau-Zener model. The model clearly shows how the 
Penning ionization process leads to the observation of the 
PAS resonances through ionization. 

Finally we apply our multi-channel model to investi- 
gate PAS of the He(23Si-23Si) system. In Ref. 1^ it has 
been shown that all observed resonances in this system 
can be attributed to the coupling of the "ground" state 
channel with four excited-state channels. We solve the 
multi-channel coupled equations for this system and ob- 
tain good agreement with the experimental results. This 
shows that the insights obtained from the simple, peda- 
gogical model can be applied to this complicated system 
of many coupled channels. Although the results we show 
are specific for the metastable helium system, they can 
easily be generalized to the other metastable rare-gas sys- 
tems and also to other systems, where ineleastic processes 
play an important role. 



II. EXPERIMENTAL SETUP 

The experimental setup consists of two main parts, the 
metastable helium source and the magneto-optical-trap 
(MOT) chamber, which are separated by three differen- 
tial pumping stages. Metastable He(2^S) and IIe(2'^Si) 
atoms are produced in a DC discharge by electron im- 
pact 7]. A Pyrex glass tube with an orifice of 1 mm 
inserted in a Teflon holder is placed in a copper assem- 
bly. The glass tube and nozzle assembly are placed in a 
cryostat. The gas flows on the outside of the glass tube 
along the cold cryostat and cold nozzle, before it expands 
into the vacuum chamber. Part of the gas is pumped out 
through the glass tube, as in the liquid nitrogen cooled 
DC-discharge source described by Kawanaka et al. @- 
However, in contrast to that source, in our design the 
source is cooled to a much lower temperature of about 
10-15 K by liquid helium. The nozzle is made out of an 
aluminum plate of 1 mm thickness to reduce heat gra- 
dients (a copper nozzle did not perform as well as an 
aluminum one) and the exit hole diameter used is 0.5 
mm. As a discharge needle we use a tungsten rod with 
a diameter of 2 mm, which is sharpened at the end and 
kept in the middle of the glass tube by ceramic spac- 
ers. The optimal distance from the nozzle plate to the 
discharge needle is found to be 7 mm, ensuring both a 



high output stability and a high yield of the source. To 
minimize heating effects the source is operated at a low 
discharge power of typically 35 mW (about 0.05 mA) 
and a low discharge pressure of 10~^ mbar. The helium 
discharge is localized inside the source, i.e., between the 
discharge needle and the nozzle plate. This is in con- 
trast to conventional metastable sources, which operate 
at higher pressures (a few mbar) with a discharge burn- 
ing outside the nozzle, i.e., between the discharge needle 
and the skimmer, to reduce quenching of the metastable 
atoms 0, ^3 ■ We observe that operation in this manner 
is only possible at higher source pressures and currents, 
and therefore at higher temperatures. However, we es- 
timated the quenching of the metastables in our design 
to be of minor importance due to our low source pres- 
sure of 10~^ mbar. Operating at minimum power and 
pressure yields a total flux of a 5x10^^ atoms/s sr. The 
mean velocity of the atoms is only 300 m/s, compared to 
2000 m/s in the case of conventional sources. Since the 
slowing distance for laser cooling is proportional to the 
square of the initial velocity, this leads to a reduction in 
slowing distance of nearly a factor 50 . 

Therefore, the slow atoms produced can be loaded 
directly in a magneto-optical trap (MOT) without the 
need of a Zeeman slower. Usually a He(2''Si) MOT is 
loaded from a discharge source cooled with liquid nitro- 
gen 1^0,01 . In that case, the atoms need to be slowed 
in a Zeeman slower down from 1000 m/s to typically 60 
m/s, before they can be captured in the MOT. Using 
a Zeeman slower makes the experiment more complex, 
since special care has to be taken to avoid reduction of 
density of the beam at the end of the slower 131 . In a 
Zeeman slower the beam is slowed down in the longitu- 
dinal direction and as it is slowed down, the transverse 
spreading of the beam becomes more dominant, thereby 
decreasing the flux in the center of the beam. This ef- 
fect limits the number of atoms that can be captured in 
the MOT. This can in principle be compensated for by 
inserting a transverse cooling section at the end of the 
slowing unit, but this will complicate the experimental 
setup further. 

To slow the atoms from the source at 300 m/s to the 
capture velocity of the MOT (« 30 m/s), we slow the 
atoms in the MOT chamber. The coils that generate the 
MOT magnetic field in our setup are large in diameter 
(about 400 mm). These large coils produce a magnetic 
field that is sufficient for Zeeman slowing and thus no ad- 
ditional Zeeman slower is needed. The atoms are slowed 
by a counterpropagating slowing laser detuned 50 MHz 
below the He(2'^Si- 2'^P2) transition with a saturation 
parameter of sq — 400. The beam still fans out in the 
slowing process, but the effect is smaller than for a liq- 
uid nitrogen cooled source for two reasons. First, due 
to the lower initial velocity of the atoms the beam does 
not spread out as much during the deceleration. Second, 
since the atoms are slowed down in the MOT chamber 
the beam spreads out mostly in the region, where the 
MOT trapping laser beams are present. This reduces 
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the losses considerably. The increased loading efficiency 
fully compensates for the low flux of the source and we 
find that the loading rate of the MOT from this source 
and the liquid nitrogen source [T]| are comparable. 

The MOT laser beam is expanded by a telescope to a 
beam with radius of 10 mm and the saturation param- 
eter of the laser light after the telescope is sq = 180. 
The laser frequency is locked by using the saturated ab- 
sorption spectroscopy 14 MHz below the He(2^Si- 2^P2) 
transition. We trap about 3x10^ atoms, as measured by 
releasing the atoms onto micro-channel plates (MCP) lo- 
cated 70 mm from the center of the MOT chamber. The 
average density of the trapped atoms is 10^ atoms/cm'^ 
and the temperature is 1.9 mK, as measured by time-of- 
flight of atoms from the MOT, using the MCP. 

An overview of the probe laser setup is shown in Fig.^ 
The probe laser beam is overlapped with the MOT cloud. 
To obtain a high probe light intensity, the probe beam 
can be focused by a lens. Additionally, a one-to-one tele- 
scope is used to overlap the focus of the probe beam with 
the MOT cloud. In order to calibrate the frequency scale, 
two beam splitters each split off 5% of the light from the 
main laser beam. One beam runs through a Fabry-Pcrot 
interferometer (FPI) with a free spectral range of « 1 
GHz and the transmitted signal is detected by a photo 
diode. This is used as a relative frequency reference. 
The other laser beam runs through a saturated absorp- 
tion spectroscopy setup, which is used to calibrate the 
offset of a frequ ency scan with respect to the He(2'^Si- 
23P2) transition [isj . 

In the experiment, pairs of cold atoms can absorb a 
photon from the probe laser beam and be excited to 
an long-range attractive potential. The probe laser is 
scanned in frequency below the D2 line and the ionization 
signal caused by Penning collisions in the MOT is studied 
as a function of the probe laser frequency. The peaks ob- 
served in the ion rate are attributed to vibrational states 
of the excited molecules. Close to the atomic D2 reso- 
nance we observe a large dip in the ionization spectrum, 
which corresponds to the fact that close to resonance 
atoms are heated due to the radiation pressure of the 
probe laser, in which case no atoms could be trapped or 
photoassociated in the MOT. 

To reduce the very high count rate of Penning ioniza- 
tion due to the MOT lasers, we periodically modulate 
the MOT laser (with a frequency 25 kHz) with a square 
wave modulation, such that during the "off" period the 
MOT laser is detuned by 500 MHz below the resonance 
frequency. This increases the signal-to-noise ratio consid- 
erably. The spectra are recorded during the trap "off" 
period and in general 20 spectra are averaged. 

We have measured photoassociation spectra for various 
probe laser intensities. In Fig. |21a number of spectra are 
shown which have been measured with increasing probe 
laser intensities. The intensities range from sq = 50-10^. 
The figure shows that, when the laser intensity increases, 
ionization peaks start to appear at increasing detunings 
and that the heights of the peaks become larger. Deeply 
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FIG. 1: Setup of the probe laser. The MOT is overlapped 
with the probe laser, which is nearly coUinear with one of 
the MOT laser beams. A small fraction of the probe laser 
light is split off for the Fabry-Perot interferometer (FPI) for 
relative frequency calibration, and for the saturation spec- 
troscopy setup for a frequency reference with respect to the 
He(2^Si- 2^P2) transition. The ions produced are detected 
with multi-channel plates (MCP). 



bound states can only be measured with high probe laser 
intensities. The bound states that lie close to the dissoci- 
ation limit are measured with low intensities since a high 
intensity probe laser would disturb the cloud of atoms 
in the frequency region of interest. We have measured 
peaks at detunings down to « 14 GHz below the D2 line. 



III. THEORETICAL MODEL 

A. Ground molecular state 

In the experiment, atoms are prepared in the He(2^Si) 
state. This state is metastable with a lifetime of 8000 s 
and thus stable on the time-scale of the experiment. We 
therefore consider it to be the "ground" state. Two atoms 
colliding in the "ground" state interact through the 
•^1]+ or potential. We only consider in our analysis 
the state based on the fact that its multiplicity is 
much higher than that of the state, and that the 
■^S^ state is repulsive. Because of the low temperature 
we also only consider s-wave scattering (/ = 0). 

In our calculations we use the long-range potential 
as calculated by Starck et al. ,1(5], given by the following 
analytical fit of their numerical calculations (atomic units 



4 



CO 

'c 

3 

SB. 

B 
fs 
cc 

c 
o 




-0.5 
A [GHz] 



CO 

B 
fo 

DC 

c 
o 



S = 10 






-4 -3 -2 
A [GHz] 



-1 



FIG. 2: Spectra measured with increasing probe laser inten- 
sity ranging from so = 50-10^ . Note, the change in detuning 
scale between the two panels. 



are used in the description of the theoretical model) 
where 



r^r^ C2nT^2n[c6{r - C7) 



n— 3 



T2n{x) - 1 



E' 2n „ 
Tn 



(1) 



(2) 



and the constants are given by ci — 0.247643, C2 
= 6.9178xl0-^ C3 = 5.851244, C4 = -0.262988, cg 
= 0.166841, C6 = 0.770951, 07 = 1.335894, = 
3.265x103, Cs = 2.075xl0^ Cio = 2.107x10^ and 
C2,i>io = C'2„-6(C2n-2/C'2„_4)^- In their paper the au- 
thors state that they use p = 5. However, we obtained 
the best fit to their data if we allowed the sum to go up 
top = 8. 



B. Excited molecular states 

The probe laser is detuned below the He(23Si-23P2) 
transition. There are many potentials that are connected 
to this asymptote, which have previously been calcu- 
lated by Mastwijk et al. and more recently by other 
groups [l7ljla|. In a joint paper by the ENS group and 
our group [g it has been shown that only four of these 
potentials are relevant for the photoassociation, indicated 
by (twice), lu and 2^ (Hund's case (c) notation). In 
Figure|3|the five potentials used in the model are shown. 
The Oj potentials are connected at short range with the 
potential, which ionizes with a probability pi of ~ 
80%. The other two states connect at short range to 
the potential, which is not ionizing. The splitting 
between these potentials is due to the fine structure cou- 
pling. The couplings between these excited-state poten- 
tials were calculated a priori and included in the model. 
It should be noted that there is no direct coupling be- 
tween the ionizing 0+ states and the 2„ state. That 
means that the 2^, channel can only couple to the ion- 
izing channels via the 1„ channel. 



C. Extended Numerov method 

The individual channels in our system without cou- 
pling arc described by the radial Schrodinger equation, 
given by 



d^ipjr) 
dr"^ 



+ k'^{r)ip{r) = 0, 



(3) 



where fc^(r) = 2pi[E — V{r)] and V{r) is the interatomic 
potential. Furthermore, E is the relative kinetic energy 
of the atoms and /i is the reduced mass of the system. To 
perform the numerical calculations presented in this arti- 
cle, we use an extended version of the Numerov method 
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FIG. 3: The ground and four excited-state potentials (solid 
lines), lu (dashed line) and 2u (dotted line). The excited- 
state potentials are shifted down by the 2^Si-2'^P2 resonance 
frequency. 

to solve the radial Schrodinger equation 'J^. In its orig- 
inal form, the Numerov equation is a discretized version 
of the one-channel Schrodinger equation. The integra- 
tion domain is divided into N steps of size h, such that 
fn = fo + nh. In that case, the Numerov equation can 
be written in a two-point form, which can be numerically 
integrated faster ^j. In the multi-channel calculations 
the wavefunction is written in matrix form, where each 
row represents an asymptotic eigenstate of the system 
(channel) and each column a linearly independent so- 
lution. The Hamiltonian contains off-diagonal elements 
that describe the couplings between these channels. For 
a model system consisting of a ground state (s) and two 
excited-state channels (pi, P2) with a r-independent ro- 
tational coupling J between the excited states, the dif- 
ference equation becomes 

F„+i = (21 + 10T„)V^„ - F„_i, (4) 

with 

T„ = - — h'' fin kp,\r.n) J . (5) 
\ J kp,^ir„)J 

Here is the Rabi coupling between the ground and 
excited states, I is the identity matrix, F„ = (I — Tn)i'n 
and fci(r„) is the wavevector of state i at r„. In this 
equation the wavefunction occurs explicitly, so inversion 
of (I — T„) is needed at every integration step. 

We start integrating outwards with some initial condi- 
tion of the components of the wavefunction of the closed 
channels. Once we pass the Condon point an excited- 
state channel becomes closed. Numerically integrating 
the wavefunction in a closed channel gives rise to two 
solutions; a physical solution, that is exponentially de- 
caying, and a non-physical solution, that is exponentially 
growing. We exploit a diagonalization procedure ^20] to 
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FIG. 4: Excited-state wavefunction. The vertical line indi- 
cates the last point where we applied the diagonalization pro- 
cedure. 

transform away the non-physical solution. The proce- 
dure works as follows: At a point r„i, where the absolute 
value of one component of F„j reaches a critical value 
that is much larger than unity (for example 10), which 
indicates that the wavefunction has grown significantly 
and therefore that the non-physical solution has become 
dominant in that channel, the inverse of Fm is calculated 
and applied to all F„ with n < m. This selects the ini- 
tial conditions such that after integration F,„ becomes 
equal to the identity matrix. With respect to before the 
inverse of F^ is applied, the norm of the wavefunction 
at Tm has decreased significantly. In fact, it is now such 
that we can conclude that the non-physical solution has 
not become dominant at r^- Since the non-physical is 
exponentially growing, this means that it does not have 
a significant amplitude at a short distance before r„i. So 
effectively we have transformed away the non-physical 
solution to select the "right" physical initial conditions. 
After continuing with the integration, we apply the dia- 
gonalization procedure again whenever the critical value 
is reached. In this way we tune the initial conditions such 
that the non-physical solution does not become dominant 
for increasing r. The result of this procedure can be seen 
in Fig. ^ which is an example for an excited-state wave- 
function. From the figure it is clear that the wavefunction 
is physical up to some small distance before the last point 
where the diagonalization procedure is applied. 

D. Accumulated phases 

For metastable helium the molecular potentials in the 
short-range domain (r <20ao) are not known accurately 
enough to use them for our integration procedure. There- 
fore we employ the accumulated phase method |2]| . In 
this method the effect of the short-range potential is 
replaced by a pre-determined phase, which represents 
the phase accumulated over the short-range region. The 
value of this accumulated phase will depend on the details 
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of the inner region of the potential, which are unknown. 
However, if the potential is sufhciently deep compared 
to the energy of the bound states, this phase is inde- 
pendent of the exact position of the bound states and 
thus one accumulated phase is sufficient to calculate all 
bound states close to the dissociation limit. As a re- 
sult, the accumulated phases can be fitted to match the 
experimental results. In the system under study some 
of the channels are ionizing in the inner region through 
Penning ionization. We take this into account by taking 
the accumulated phase to be a complex quantity in this 
case, where the real part is the actual phase shift of the 
wavefunction and the imaginary part takes into account 
the ionization. 



F. Extracting cross sections 

The key to extracting cross sections from the calcu- 
lated wavefunctions is the phase shift 6i in the open chan- 
nel (ground state). This phase shift describes the asymp- 
totic effect of the collision between the atoms. Asymp- 
totically the component of the wavefunction in the open 
channel can be written as 

Rtik,r) = Ae{k) {je{kr) - tan Si,{k) m(kr)) , (8) 

where i is the partial wave, and rji are the spheri- 
cal Bessel and Neumann functions, respectively, and the 
Ai(k) are prefactors. As mentioned earlier, in our numer- 
ical calculations we only consider ^ = Q. In the asymp- 
totic region the solution for the open channel reduces to 



E. Integration 



ipsir) = A{ks) cos So sin(fcsr + (5o), 



(9) 



For the numerical integration of the wavefunction, the 
overall error in the integration is important. Blatt p3 | 
derived that for the Numerov method the relative error 
per step £i can be approximated by 



Error in ^> 



n+l 



72, 
10 



T^oih^'k'' r^K'lE-Vir)]^ 



(6) 

However, as shown by Sloan |23j | the error ejv after many 
steps becomes proportional to h^. We find that 



£n --h'^\E-V{r) 



(7) 



This shows that it will be advantageous to have our step- 
size dependent on — y(r)| with a smaller stepsize for 
the region where V{r) is still very deep. This is the basis 
for the Fourier grid method However, in our case, 

where we are dealing with a few coupled equations, we 
restrict ourselves to changing the stepsize hi = h once to 
a new stepsize /i2 = fcft-i, where fc is a positive integer. 

In practice we use 100.000 steps for the one-channel 
calculations to find the bound states in the potentials. 
The integration range runs from 20 to 900ao. Half the 
number of steps are used for the first 80ao, whereas the 
other half of the steps are spread out over the remaining 
SOOao. For the multi-channel calculations we integrate 
from 20 to 2000ao using 10.000 points. Half these points 
are equally spread between 20 and 60ao, and the other 
half is spread out over the rest of the integration range. 
Depending on the detuning, the Condon points for these 
four potentials generally lie between 120 and 240ao in our 
calculations. Also depending on the detuning, we need to 
apply the diagonalization procedure somewhere between 
10 and 100 times. The total flux is then conserved over 
the integration range to an accuracy of more than 98% 
for the final calculations presented in this paper. This is 
an indication of the accuracy of the procedure. 



where ks is the wavenumbcr for the s-channcl in the 
asymptotic region, where it is no longer dependent on r. 
The logarithmic derivative of this last equation is given 
by 



7(r) 



1 di}}s{T') fcs cos(fcsr) — fcs tan(5o sin(fcsr) 



ips{r) dr 



sin(fcsr) — tanJo cos{ksr) 



We can solve this for tan 5q to find |25| 

kg cos(fcsr) — 7(r) sa\{ksr) 



ks sin(A;sr) -|- 7(r) cos(fcsT') 



(10) 



(11) 



The logarithmic derivative in the Numerov method is 
given by 

7„ ^ ^ = i f A„+ii?„ - ^\ (1 - r„), (12) 

where An = {\ — T'„)/(l — T„). We now use these two 
equations to determine the phase shift 5o in the open 
channel. The point where we determine the phase shift 
should lie a short distance before the point where we last 
applied the diagonalization procedure to ensure that the 
non-physical solution does not contribute to the results. 
Since we consider ionization in the photoassociation pro- 
cess, the atoms can scatter either elastically or inelasti- 
cally, which in our case is ionization. Using the complex 
phase shift (5o = ^o,i?, + *<^o,/ we denote the elastic cross 
section UeUs as (25.J 



47r 
fc? 



sin^ (5o,i?, 



and the inelastic cross section dincias as 



^inclas 



fc? 



[1 



-45o,. 



(13) 



(14) 



The last expression will be used for the identification of 
the ionization process in our experiment. 
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FIG. 5: The uncoupled "S+S" and "S+P" molecular states 
(dashed lines) cross at the Condon point Rc leading to an 
avoided crossing of the coupled states (solid lines). 



IV. MODEL CALCULATIONS 
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FIG. 6: Breit-Wigner profile in the elastic cross section (log- 
arithmic scale) as a function of the detuning 5 in the two- 
channel case with C3 = -10, E = and H = 66 MHz. 
The crosses are numerical data and the line is a scaled Breit- 
Wigner profile with Er = -28.4374 GHz, F = 998.8 kHz and 
cj} = -3.721x10"^ rad. 



A. Introduction 

In this section we examine in detail some simple model 
systems that illustrate the physics of the ionization pro- 
cess in the metastable helium experiment. A theoretical 
introduction is given before we describe the model sys- 
tems. The model systems consist of a flat ground-state 
potential and one or two excited-state potentials. In both 
cases one excited-state potential is ionizing. Where appli- 
cable, we compare our numerical results to Landau-Zener 
type calculations. 



B. Avoided crossing 

The potentials in the helium system were calculated 
using the Born-Oppenheimer approximation and includ- 
ing the effects of rotational coupling between the different 
states 0, . The long range potential in the "ground" 
("S-l-S") state is a van der Waals l/r^ potential and in 
the range of interest to us it can be considered flat. The 
excited-state ("S-l-P") potential has a I/r"^ behaviour at 
long range. For simplicity, we add the energy of one 
photon to the ground-state potential. For a negative de- 
tuning 5 the shifted potentials cross at the Condon point 
Rc^ as shown in Fig. O Due to the coupling of the hght 
field such a crossing becomes avoided and the potentials 
in the vicinity of the crossing are modified. 



C. Two channels without ionization 

In this model we neglect for simplicity the fine struc- 
ture interaction in the 2^V2 state. Traditionally such 
a two-channel system can be approximately described 
analytically using the Landau-Zener method. Using the 



approximations discussed above, we describe the poten- 
tial for the "S-fP" state as Vs+p{r) = 6 + Cz/r'^ with 
C3 the dispersion coefficient for the dipole-dipole in- 
teraction. The Condon point is in this case given by 
i?c - ("C3/<5)i/^ 

An example of a photoassociation resonance is shown 
in Fig. IHl where we show the elastic cross section for the 
non-ionizing case. The elastic cross section consists of 
two parts. The first part is an almost constant back- 
ground, which arises from the background scattering in 
the ground-state potential. The second contribution is a 
resonance when the energy of the incoming atom matches 
with a bound state in the excited-state potential. The 
two contributions interfere coherently to form a Breit- 
Wigner profile 25] given by 



TT 



r/2 



{5-E-E.r)+iT/2 



(15) 



with Er the position of the resonance, F the linewidth 
of the resonance and (p the phase difference between the 
direct scattering (the background) and the resonance. 

To take a closer look at what goes on at a photoassoci- 
ation resonance, we consider Fig. [7Ji. On the right-hand 
side of this figure the bound states in part of the energy 
spectrum of the C3 = —10 potential for an accumulated 
phase of rad at 20ao are indicated. We assign numbers 
to these bound vibrational states by counting the nodes 
of the wavefunction. However, since we start to integrate 
at 20ao, the number of nodes in the inner region is un- 
known. As a result a state numbered n is not the n*'' 
vibrational state of the potential, but rather a state with 
number n -|- m, where m is an offset. On the left-hand 
side of Fig.Et we show a plot of the maximum amplitude 
of the excited-state wavefunction in the two-channel cal- 
culation. It is clear from the graph that the wavefunction 
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maximizes, whenever the energy corresponds to a bound 
state. 

Some wavefunctions are drawn in Figs. 03 and Et. 
From Fig. [7|d it is clear that the effective couphng for 
the three wavefunctions around the resonance is much 
stronger than for the other two, as their phases are al- 
tered much more dramatically. From Fig. |3; we can see 
that below the resonance (the dot-dashed line) we indeed 
have one node less than at resonance or above it. Also 
the artificial, exponential increase of the wavefunctions 
due to the renormalization procedure near the end of the 
integration range is visible. 

In Fig |S1 we have plotted the resonance position Er , 
the width F and the phase difference as a function of 
the accumulated phase (j)s of the "S-l-S" potential at 20ao, 
which we obtained by fitting the calculated cross sections 
using Eq. E| If we change the phase of the "S-f S" po- 
tential, the Franck-Condon overlap between the "S-l-S" 
and "S-l-P" wavefunction is changed (see Fig.^. If the 
"S-|-S" wavefunction has a node at the Condon point, 
the overlap between the two wavefunctions is small. This 
leads in Fig.lHb to a very narrow resonance for (j>s = O.Ott, 
whereas for 0s = 0.457r the overlap is maximal and thus 
the width is the largest. The coupling not only leads to 
a broadening of the resonance profile, it also leads to a 
shift of the resonance position (see Fig.|S^). This shift is 
called the light shift, since it is induced by the coupling 
of two potentials by the laser light. The sign of the light 
shift depends on the fact, whether the overlap is stronger 
for distances smaller or larger than the Condon radius 
Rc- For smaller distances the "S-l-P" potential energy 
is lower than the "S-l-S" potential energy and thus the 
couphng between the two potentials pushes the "S-l-P" 
potential down in energy, which leads to a negative light 
shift. Conversely, if the overlap is larger at distances 
larger than Rc, this leads to a positive light shift. Sur- 
prisingly, if the overlap between the two wavefunctions is 
maximal and the coupling is the strongest {(f>s = 0.457r), 
the light shift becomes zero. This is caused by the fact 
that the "S-f S" wavefunction peaks at the Condon point 
and the "contributions" to the light shift for smaller and 
larger distances cancel. Finally, we see in FiglS); that the 
phase difference is directly proportional to the accu- 
mulated phase (j)s , which is to be expected in this simple 
model. 



D. Two channels with ionization 

We extend our model by adding an ionization process 
in the inner region of the excited-state potential and con- 
sidering also the inelastic cross section. An example of a 
resonance in the inelastic cross section is shown in Fig. 1101 
Exactly as the elastic cross section, it contains resonances 
from bound states in the excited-state potential, which 
take the form of Lorentz profiles due to the absence of a 
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FIG. 7: Right-hand side of Fig. (a) is an energy diagram. On 
the left-hand side a logarithmic plot of the maximum value 
the excited-state wavefunction reaches. Figure (b) is a plot 
of some ground-state wavefunctions around the n — 21 reso- 
nance. The solid line is for S slightly above the resonance, the 
dashed line is calculated at the resonance (—28.4374 GHz), 
the hashed line at the zero-crossing of the resonance, the dot- 
ted line at the lowest point of the resonance, and the dot- 
dashed line is for well below the resonance. Figure (c) is the 
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FIG. 8: Resonance properties as a function of the accumu- 
lated phase (f>s of the "S-|-S" channel at 20ao: (a) Position 
of the resonance Er, (b) width of the resonance F and (c) 
phase difference (j) between the direct scattering and the reso- 
nance. For small couplings the position of the resonance is at 
—28.43543 GHz. Due to the changing Franck-Condon overlap 
between the ground and excited state, both the shift and the 
width of the resonance change. 
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FIG. 9: Wavef unctions of the "S-l-P" channel (solid line) and 
"S-l-S" wavefunction (dashed, dashed-dotted) for two differ- 
ent accumulated phases of the "S-l-S" channel. For the first 
value of the accumulated phase (dashed line) the wavefunc- 
tion has a node at the Condon point {Rc ~ 160ao) and the 
overlap between the two channels is small, whereas for the 
second value of the accumulated phase (dashed-dotted line) 
the overlap is ma^ximal. 
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FIG. 10: Lorentz profile on a logarithmic scale in the inelastic 
cross section as a function of 5. The crosses are numerical data 
and the line is a Lorentz profile with Er = —28.47106 GHz, 
r = 4.770 kHz and A = 1.583x10^ ao^ 



direct ionization signal, given by 



'^inela 



(16) 



where A is the height of the Lorentz peak, F is its width 
and Er is its position. 

For comparison, we develop a semi-classical picture, in 
which a metastable helium atom comes in on the "S-l-S" 
potential. At the Condon point, wc denote with pc the 
probability that the atom follows the adiabatic route to 
the "S-l-P" potential. The probability that the atom ion- 
izes in the inner region of the "S-l-P" potential is des- 
ignated with Pi. From the inner regions of the two po- 
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FIG. 11: Flux diagram for the two-channel problem with ion- 
ization. At each point (small circles) the system has a certain 
probability to make a transition at the crossing. In the end 
the system either scatters elastically (open circles) or inelas- 
tically (filled circles). 



tentials, and from the outer region of "S-l-P", the atom, 
if not ionized, is assumed to reflect back to the Condon 
point. In Fig. 1111 a flux diagram is drawn for this semi- 
classical situation. It contains three simple loops (i.e., 
loops without external branches) that can be integrated 
out. The total probability Pi for the atom to ionize is 
then found to be 

p ^ l + (l-pc)V(l-pg) ^ 2p,p, 

(17) 

In the middle part of this equation the three geometric 
series for the three loops can still be recognized. We can 
calculate the crossing probability Pc using the Landau- 
Zener formula 

Pc^l-e-^^, (18) 
ZaVc 

where n is the Rabi frequency between the ground and 
excited state, a is the difference in derivatives of the po- 
tentials at the Condon point, and Vc is the velocity of 
the incoming atom at the Condon point. For a kinetic 
energy E of 1.9 mK, a Rabi frequency of 66 kHz and a 
detuning of 6 ^ -28.4 GHz, we find pc = 8.84x10^1". 

One should note that in such a semi-classical treatment 
of the problem we can no longer treat the dependence in 
6 and only the total ionization probability is calculated. 
Quantum mechanically, however, we calculate cross sec- 
tions that are a function of 6. As a measure of the total 
ionization probability, we have therefore taken the total 
area underneath one peak of the inelastic cross section. 
Although the integrated inelastic cross section obtained 
in this way is proportional to the total ionization proba- 
bility, we lack a method to map the quantum-mechanical 
results onto the semi-classical picture. Therefore we can 
only compare the two methods with respect to their qual- 
itative behavior as a function of the ionization probability 
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FIG. 12: Comparison between semi-classical and numerically 
calculated values for the total ionization probability Pj as a 
function of the ionization probability pi in the excited-state 
potential. The crosses are the scaled numerical data, the solid 
line is Eg. 1171 with a fitted value for pc (non-scaled), and the 
dashed line is the same equation with the Landau-Zener value 
for Pc and scaled. The triangles are numerically calculated 
widths r of the resonances, and the squares are their heights, 
both in arbitrary units. 

Pi, but not with respect to the absolute total ionization 
probability Pj they give rise to. 

In Fig. ^] we plotted the results of both calculations as 
a function of the probability pi. We normalized for large 
Pi the total amplitude of the numerical calculation to the 
total ionization probability calculated semi-classically. If 
we allow Pc to be an adjustable parameter, the best fit 
is for Pc — 4.59x10^^°, which is about half the Landau- 
Zener value. Note, that the Landau-Zener formula de- 
scribes a single crossing, whereas in the current situation 
the atoms vibrates back and forth over the crossing lead- 
ing to interference effects. From the figure it can be seen 
that the widths of the resonances are increasing as a func- 
tion of Pi, but their heights reach a maximum around the 
point where P/ reaches its asymptote. 



E. Three channels with ionization 

Next we consider the case where we have two excited- 
state potentials. We take both excited potentials to be- 
long to the same "S-fP" asymptote, but have slightly 
different C3 values. These potentials will be denoted 
"S-HPl" (C3 = -9) and "S-I-P2" (C3 = -10). Only the 
"S-l-P 1" potential leads to Penning ionization at short in- 
ternuclear distance. The potentials are coupled by a ro- 
tational coupling with a probability pj , which acts over a 
large distance. This system is depicted in Fig. E| Since 
the effective coupling between two potentials is propor- 
tional to the wavefunction overlap and the wavefunctions 
of the closed channels have a maximum close to their 
Condon points, we expect the rotational coupling to be 
most effective before the first Condon point. For the 
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FIG. 13: The model system with two excited-state potentials. 
Only the "S+Pl" potential is ionizing with probability pi. 
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semi-classical treatment of the total ionization probabil- 
ity P/, we assume that the coupling only takes place at 
this point, and that the atom can make a transition from 
one "S-l-P" potential to another both on the way in and 
on the way out. Furthermore, in the Landau-Zener model 
the transfer probability depends quadratically on the off- 
diagonal term in the Hamiltonian. 

In Fig.^lwe consider the resonance at 6 = —28.4 GHz, 
where we assume that the two states are coupled to the 
ground state with the same Rabi frequency. Since the 
Rabi frequency is large, the resonances from the "S-I-Pl" 
potential are so wide that they overlap to a large extent. 
In contrast, bound states in the non-ionizing "S-fP2" 
potential will give rise to much sharper peaks in the in- 
elastic cross section, because only the part of the wave- 
function coupled to the ionizing "S-I-Pl" potential can 
lead to ionization. The inelastic cross section thus con- 
tains a background part caused by direct ionization in the 
"S-I-Pl" channel, which interferes with resonances from 
the bound states in the "S-I-P2" potential. This leads to 
a Breit-Wigner profile with a non-constant background. 
The result of this interference can be quite unpredictable 
and instead of fitting the profiles we choose to numer- 
ically integrate the area underneath the resonance as a 
measure of the ionization probability. 

To solve the semi-classical problem, we draw a flux dia- 
gram containing all the possible paths of the system. For 
every intersection in this diagram we write down an equa- 
tion, and combine these equations into a set of coupled 
linear equations that is solved analytically. The solution 
for the total probability Pi is given by 

P _ 2k ((1 - Pi)piP2 + xjpj - p])) 

rj , (iyaj 

Px+Py+ Pz 

Px = Pi (P2[2 - 2pi(l -pi) -pi] +pi -plPi) , (19b) 

Py = {2x+ (1 -p2 -pi[4- 6j92 -Pi(3 - 7p2)])pi)pj, 

(19c) 

p, = - {2x + (1 - 2pi)[l - 2p2-pi{l - 3p2)]p^)pl 

(19d) 

x^P2 + {l-pi)pi{l-5p2). (19e) 
For the parameters used in Fig. 1141 the Landau-Zener 



FIG. 14: Resonances in case of a single kinetic energy (dashed 
line) and a thermal distribution (solid line). The Rabi fre- 
quency is f2 = 6.6 MHz for both couplings to the excited state, 
an ionization probability pi = 0.80 in "S-I-Pl" and an energy 
of E = 2.0 mK. For the thermal averaging the temperature is 
E/kg. 



probabilities are given by pi = 8.84x10^^ and p2 = 
8.54xl0~^, where the small difference stems from the 
fact that the "S-I-Pl" potential is less steep at its Con- 
don point than the "S-I-P2" potential. 

In Fig. 1151 we compare the semi-classical and numeri- 
cal calculations as a function of the rotational coupling 
Pj. The value of Pj for large pj is determined by pi, 
whereas for small pj it is determined by p2. To get good 
agreement between the numerical analysis and the semi- 
classical description we have adjusted the parameters pi 
and p2 and find pi = 9.12x10"'' and p2 = 3.60x10"''. 
Note that the relatively large difference between these 
two parameters, which are nearly identical in the Landau- 
Zener description (see above), reflects the fact that for 
these two states the Franck-Condon overlap with the 
"ground" state at the Condon point can be very differ- 
ent and thus that the behavior is dominated by Franck- 
Condon factors. Although the factors pi and p2 calcu- 
lated this way do not each have physical significance, 
their ratio Pi/p2 has. Numerically, we calculate this ra- 
tio to be 2.53, whereas Landau-Zener description yields 
1.04. Clearly, we have to resort to quantum mechani- 
cal methods to solve this model, since the Landau-Zener 
formula assumes isolated avoided crossings, whereas we 
have two crossings lying close to each other. From Fig. 1151 
it is clear that with the fitted parameters pi, p2 and pj 
the agreement between our numerical calculations and 
the semiclassical result is very good. 



F. Thermal distribution 

In our experiment, the interacting atoms do not all 
have the same relative energy, but instead have a thermal 
distribution determined by the temperature of the sam- 
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FIG. 15: Numerical data for the three-channel system as a 
function of the rotational probability pj (crosses) and the 
semiclassical formula of Ea. ll9l fitted to the data (line), where 
we fit the semiclassical parameters pi, p2 and the coupling 
strength between pj and the square of the off-diagonal matrix 
element. 



pie. We investigate the effect of the thermal distribution 
by integrating the profiles over a Maxwell-Boltzmann 
thermal distribution. The effect of the thermal averaging 
depends on the width of a single resonance compared to 
the Doppler shifts. In Fig.^lwe show a plot of the reso- 
nance at i5 = —28.4 GHz in the inelastic cross section for 
a relative energy E = fc^T and for a thermal distribu- 
tion. In the latter the dip has largely disappeared. This 
is as expected, as the Doppler width and the resonance 
width are about the same size here (^6x10^^ E/j). 



V. COMPARISON BETWEEN THEORY AND 
EXPERIMENT 

Now we return to the numerical treatment of the exper- 
iment. As is clear from the previous section, the location 
of resonances in the inelastic cross section depends on 
the accumulated phases in the non- ionizing (1„ and 2^) 
channels. In order to determine the accumulated phases, 
we first have to identify the resonances belonging to one 
particular potential. In a combined effort between our 
group and the group at ENS in Paris P| this has been 
done for this system. This effort is based on the fact 
that all the resonances in one potential have (at least 
approximately) the same accumulated phase. So by cal- 
culating the accumulated phases needed to reproduce all 
resonances, one can group the resonances together. 

As mentioned before, we are left with two sets of res- 
onances; one for the 1„ potential and one for the 2^ po- 
tential. Table is a list of the resonances in the two 
channels. For both of these potentials we calculate the 
bound states as a function of the accumulated phase in 
a simple one-channel calculation. Selecting the accumu- 
lated phase that correctly reproduces one resonance of 
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11.10 


13.57 


7.01 


8.94 


4.26 


5.64 


2.42 


3.38 


1.21 


1.88 




0.98 




0.46 




0.19 



TABLE I: Part of the set of resonances (J — 3) in the two 
non-ionizing potentials being considered. 

that potential indeed reproduces all the other resonances 
of that channel at the correct positions. The phases cal- 
culated using this procedure are O.SGStt for the 1„ po- 
tential and 0.4297r for the 2„ potential. These turn out 
not to be the entirely correct phases to use in the multi- 
channel calculations, especially so for the 2„ potential. In 
general, in calculations performed for the simple models 
of the previous section, accumulated phases for one and 
multi-channel calculations matched as expected. That 
this is not the case here is due to the strong coupling 
between the excited-state potentials, which causes the 
resonances to shift. This explanation is supported by 
model calculations for the three-channel model, in which 
the resonances are observed to shift for strong rotational 
coupling between the excited-state channels. Through 
fitting with experiment we find the accumulated phases 
to be 0.3637r for the 1„ potential and 0.4457r for the 2„ 
potential. This means that the rotational coupling has 
an effect on the exact location of the resonances. The 
accumulated phases of the two ionizing excited-state po- 
tentials have no influence on the spectrum and we simply 
take these phases equal to zero. 

A. Accumulated phase ground state 

In the previous section we discussed the selection of the 
appropriate accumulated phases for the excited states. 
However, also in the ground-state channel we need to 
select the correct accumulated phase. This phase has 
a very large influence on the s-wave scattering length 
as this scattering length depends directly on the final 
phase shift of the wave function. Recently Moal et al. 
measured it very accurately to 141.96 ± 0.09 ao- We 
choose the accumulated phase for the ground state to 
match this value using |25j 

a = -lim-;^, (20) 
fc^o k 

where a is the s-wave scattering length. We use a low 
energy of 19 /iK and a low laser intensity of Sq = 10. 
That way we find the accumulated phase to equal 0.9397r. 
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FIG. 16: Figure (a) is a numerical spectrum for T — 1.9 mK 
and So — 1x10*. Figure (b) is an experimental spectrum for 
T = 1.9 mK and so a few times 10*. 



B. Spectrum 

To obtain a numerical spectrum that we can compare 
with experimental results, we use the accumulated phases 
from the previous sections and we integrate over a ther- 
mal distribution as explained in Sec. IIV Fl Figure 116b 
is calculated for a temperature of 1.9 mK and a laser 
saturation parameter of so = 1x10*. The former is in 
accordance with experiment, the latter is low in com- 
parison with experiment. A comparison with the exper- 
imental spectrum that is reproduced in Fig. 116b shows 
that we successfully reproduce the resonances at the cor- 
rect positions. Both spectra contain resonances that have 
widths in between roughly 100 and 200 MHz. However, 
the thermal effects do not cause enough smearing of the 
resonances to make the resonances symmetric in our cal- 
culations, whereas they look completely symmetric in the 
experimental data. The theoretical spectrum contains 
also a much higher peak-to-background ratio than the 
experimental data. Note that also the so-called optical 
collisions, i.e., the increase of the background level as we 
approach zero detuning, is present in our calculations. 
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FIG. 17: Two resonances at different temperatures. Solid line 
is 1.9 mK, dashed line is 2.3 mK and dotted line is 2.7 mK. 



This is because for small detunings many bound states 
are available in the excited-state channels. Therefore a 
very large number of resonances of those channels con- 
tribute to an increased ionization rate. Note, that for 
small detunings our assumption that we only have s-wave 
scattering breaks down. 



C. Line widths 

To study in more detail the effects that determine the 
line widths of the resonances in the numerical spectrum, 
we perform the same calculation as in the previous para- 
graph, but now for different temperatures. The results 
can be found in Fig. El The figure shows that a higher 
temperature will lead to a slightly narrower resonance. 
This is counterintuitive, since one would expect that a 
higher temperature leads to a larger Doppler broadening. 
This is also observed for the model systems of Sec. IIVI 
In this case the width of a resonance seems to be deter- 
mined by other effects. This is in accordance with what 
is observed in the experiment. 

One possible suggestion is that the line widths are de- 
termined by the last oscillation of the bound state at 
the Condon point. This final part of the wavefunction is 
much larger than the remainder of the wavefunction, and 
thus could therefore dominate the Franck-Condon factor. 
Its projection on the potential gives us the distance over 
which photoassociation can occur, and thereby also the 
range of detunings that contribute to the photoassocia- 
tion process. In Fig. llSl we show that at the temperature 
used in the experiment the contribution to the Franck- 
Condon factor of the other oscillations is significant. Fur- 
thermore the last oscillation covers a range of the order 
of 50ao, leading to a detuning range of 10 GHz. This is 
much wider than the resonance widths observed in the 
experiment, and we conclude that this last oscillation is 
not the limiting factor for the resonances studied. 
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FIG. 18: The wavefunction of the bound state at 5 = —8.94 
GHz depicted in its 2u potential. 

VI. CONCLUSIONS 

In this paper we have studied the role of Penning ion- 
ization for the observation of photoassociation in the 



He2(2'^Si-2^Si) system. We have shown that the reso- 
nances in experimental ionization rates are due to bound 
states in the long-range potentials near the 2'^Si-2'^P2 
asymptote that are reached through photoassociation. 
We can describe this process with the multi-channel 
model discussed in Sec where the coupling between 
bound, ionizing channels and bound, non-ionizing chan- 
nels leads to well-resolved peaks in the ionization spec- 
trum. Although the line shape is not completely un- 
derstood from the model, the location of the lines is in 
good agreement with the predictions of the model. The 
study shows that photoassociation spectroscopy is a pow- 
erful technique at ultracold temperatures even in systems 
where inelastic processes play a major role. 
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